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Physics of wave propagation and wave excitation in black hole spacetimes has been 
investigated as a tool to comprehend properties of black hole spacetimes. Especially, 
the quasi-normal mode of black holes is directly related to the structure of black 
hole spacetimes and is expected to be confirmed observationally in near future using 
gravitational waves. In the high frequency eikonal limit, the quasi-normal mode is 
connected to the wave scattering in the vicinity of the peak of the effective potential 
and thus related to the existence of the unstable circular orbit for null rays around the 
black hole. The geometric interpretation of the quasi-normal frequency is as follows; its 
real part corresponds to the period for a massless particle or null rays in the unstable 
circular orbit and its imaginary part represents the decay rate from the unstable circular 
orbit [DEIEIII]. 

On the other hand, also related to the existence of the unstable circular orbit, 
properties of “shadows” of black holes ISIEIIZ! have been studied in detail recently. The 
main motivation of this subject is direct verification of black hole spacetimes using VLBI 
radio telescopes. The shadow of the black hole is dehned as the region of absorption on 
the observer’s screen and its rim corresponds to the unstable circular orbit of null rays 
around the black hole projected on the observer’s screen. Its shape can be drawn by 
integrating null geodesics in the black hole spacetime (ray tracing). Recently, we have 
investigated the imaging of black holes using waves 019 ]; instead of applying the ray 
tracing method, the image of the black hole is obtained using the Fourier transformation 
of the scattering wave at the observer. In this method, wave effects such as interference 
can be included in images of black hole shadows. 

Wave optical treatment of the gravitational lensing has already done and 
detectability of the interference effect is discussed for the weak gravitational lensing [TOl 
HB na US). For ordinal astrophysical sources, we cannot expect to observe the spatial 
interference fringe pattern because the coherent time of the source is too short compared 
to the path difference in the gravitational lensing system. However, even for such a case, 
the interference effect can be observable in the frequency domain as the oscillation of 
the power spectrum. The period of the oscillation in the power spectrum is determined 
by the mass of the gravitational source which gives rise to the gravitational lensing na. 
For the gravitational lensing by the black hole, a new interference effect associated with 
the existence of the unstable photon orbit is expected; null rays can go around the black 
hole arbitrary number of times (orbiting) and direct rays and these winding rays can 
interfere and the beat appears in the power spectrum. We expect that the period of 
the beat is related to the mass of the black hole and the radius of the unstable circular 
orbit. 

In this paper, we consider wave optics in the Schwarzschild spacetime using a 
massless scalar field to investigate wave effects peculiar to black hole spacetimes. The 
massless scalar field is adopted as the benchmark treatment for the wave scattering 
problem by black holes and we do not consider polarization degrees of freedom which 
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is necessary for the electro magnetic waves and gravitational waves. We follow the 
standard treatment of the wave scattering problem by black holes [I51[I61[I9] and derive 
the asymptotic form of the scattered wave at sufficiently far from the black hole in the 
high frequency eikonal limit. The wave is represented by the sum over the partial waves 
which include the WKB phase shift. Comparing to the standard treatment of the wave 
scattering theory, our formulation retains the next to leading order contribution in the 
phase of the scattering wave at large r, and this rehnement makes the sum over partial 
waves converge even for the scattering caused by the long range 1/r potential. 

We shows that the gravitational lens equation in the black hole spacetime can 
be derived by taking the short wavelength limit of the scattering wave. This lens 
equation is an extension of the relation between the scattering angle and the derivative 
of the WKB phase shift which is derived by the stationary phase method [201 El] ■ To 
evaluate the sum over partial waves, we apply the complex angular momentum (CAM) 
method [221 |23l [m Hg ESj. The scattering wave is decomposed to the direct part and 
the winding part. The winding part can be evaluated using Regge poles which are 
associated with the existence of the unstable circular orbit around the black hole. 

The outline of this paper is as follows. In Sec. II, we introduce the WKB Green 
function and derive the lens equation from the Green function by the stationary phase 
method. In Sec. Ill, we shortly review the S-matrix and the Regge poles. We evaluate 
the sum over partial waves and obtain the asymptotic form of the scattering wave in 
Sec. IV. As applications, we consider images of the black hole and the interference effect 
in the power spectrum in Sec. V. Sec. VI is devoted to summary. We use units in which 
c = h = G = 1 throughout this paper. 

2. Wave scattering by a black hole 

We investigate the wave scattering problem by the black hole in the high frequency 
eikonal limit. The high frequency condition is Meo 3> 1 where M is the mass of a black 
hole. Thus this condition states that the wave length is sufficiently shorter than the 
size of the black hole. The eikonal limit is the angular quantum number i of the wave 
is larger than unity £ 3> 1. Under these conditions, we follow the standard approach of 
the wave scattering theory; decompose scattered waves to the sum over partial waves 
and introduce the phase shift which represents the properties of the scatterer. 

2.1. Wave equation and the WKB phase shift 

Let us consider the scalar wave equation = 0 in the Schwarzschild spacetime. The 
metric is 



( 1 ) 
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As the geometry has the spherical symmetry, for the wave conhguration with axial 
symmetry about the 2 ;-axis, the scalar held <h is separable as 

^ = e-^‘^^Ri{r)Pi{cose). ( 2 ) 

The radial wave function Ri and the angular wave function satisfy 


d ( . dRi\ ( 


— A 
dr \ dr 


+ 


A 


-i{i + l)]R, = 0, 


1 d 




+ ^)Pi — 0, 


( 3 ) 

( 4 ) 


sin 9 dd 

where A = r"^ — 2Mr. We consider the WKB solution of the radial equation with the 
assumption of the high frequency eikonal limit Mu 3> 1, f' ^ 1. By introducing a new 
radial function 'ipi = rP^ and the tortoise coordinate = f dr^ = r + 2M ln(r/(2M) — 
1), the radial equation becomes 
d^ijji 


drl 

Q = oj'^ - 


+ QiJi = 0 , 


( 5 ) 


+1) ^ ‘m 


1 - 


2M 


The lowest order WKB solution of this equation is 


= n-V4 




exp 


dr^ 


oc e 


iSf 


( 6 ) 


For the angular wave function, which is the Legendre function in the present case, it 
also can be written in the WKB form Pi oc For £, Mu 3> 1, the phase functions Sr 
and Sg satisfy the following Hamilton-Jacobi equations 

2 1 _ 

( 7 ) 


if-) 


dSe 

dO 


= 0 . 


( 8 ) 


These equations are the same as the Hamilton-Jacobi equations for null rays in the 
Schwarzschild spacetime. The solutions of these equations are 


Srir) = 


dr 


Vn 


5^(0) = / d9L = L9, 


(9) 


where P{r) = u‘^r^ — L^A and we have replaced i{£ -|- 1) by L'^ = {i + 1/2)^, which 
improves the accuracy of the WKB approximation; this replacement ensures the phase 
shift introduced in (16) zero for the flat limit M = 0. Following the standard Hamilton- 
Jacobi theory, by differentiating the function Sr + Sg with respect to the angular 
momentum L, the equation for the trajectory of null rays is derived as 
'"'f dr 


-L 


In 


Vn 


= -{ef-e,). 


( 10 ) 


From this equation, the deflection angle for null rays is obtained by taking rj,rj —)■ oo: 

’ dr 


0 = TT — {9f — 6i) = TT — 2L 


iro 




( 11 ) 
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where tq is the radial turning point determined by TZ{ro) = 0. 

We rewrite the WKB radial function using the phase shift. For this purpose, we 
consider the solution of the radial equation ([^ with the following boundary condition: 


U-y 


Aoute*"’'‘+Ane- 
The phase shift is defined by 

2i5i _ / \r^out 


r* —)■ —oo, 

r* —)■ +CX). 




A, 


( 12 ) 


(13) 


Then the wave at r —)■ +cx) is expressed as 

■ f r 

oc sin ( (ur* + 0 £- — 


(14) 


For 6g = 0, Uin corresponds to the spherical wave in the flat spacetime jiiur) ~ 
sin(ci;r — 'nljT). Within the WKB approximation, the analytic continuation around 
the turning point yields the radial function for large r as 




sm 


dr* ^ 


Lort 


= sm 


^ 7r£ 

wr* + di- — 


d/- ((?■/= - u) 


(16) 


where is the turning point determined by Q{rt) = 0. Hence the WKB phase shift in 
the Schwarzschild spacetime is obtained as [inj 

-2 ^ / 1' 


S£= I d/— -u) - uTt* + ^ (^ + 2 


'rt 


(16) 


For small i less than a some critical value ic, there exists no real solution for Q{r) =0 


and (16) has no meaning. However, by performing a suitable analytic extension of the 


formula, it can be shown that the phase shift acquires an imaginary part which represents 
the absorption of waves by the black hole. In the eikonal limit £ ^ 1, Q{r) cs TZ{r)/r^ 
and the phase shift is 





r \ TT 

A - ‘^>■0. + 2 


1 

+ 2 


(17) 


and by differentiating with respect to the angular momentum £, 

dSi T£ f „ 1\ dr . . 

le’^2~W2)j„ yi' 

Comparing to Eq. , we thus obtain the well known relation between the classical 
deflection angle and the WKB phase shift 


( 19 ) 
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For large r, the integral in the phase of the radial wave fnnction (15) can be approximated 
to be 


dr { 1 + — u) ^ I dr 


(£+ 1 / 2)2 (£+ 1 / 2 )' 


2ur‘^ 


2ujr 


,( 20 ) 


where we assnme that the fnnction jg go £ mnst satisfies \/2ur > £+1/2 to make 
this approximation have meaning. This imposes the maximal npper limit £max ~ tor 
for a fixed valne of ur when we take the partial wave snm. 

In the standard formnlation of the wave scattering problem, this term is treated as 
zero because it vanishes as r —)■ cx) and does not contribute to the scattering amplitude 
and the differential cross section. However, as we will see, this term is necessary 
to reproduce the lens equation from the wave function in the eikonal limit and also 
necessary for convergence of the sum over partial waves. 

After all, for sufficiently large r, we obtain the following form of the radial function 
in the eikonal limit 

with i<V2ujr. (21) 


Uin[r) oc sm 


(£+ 1 / 2)2 

tor^ -\ --- V 5^ —— 

2Lor 2 


2.2. Green function 

We aim to obtain the scattered wave by a black hole; the scalar wave is emitted from 
a point source and scattered by the Schwarzschild black hole (see Fig. 1). For this 
purpose, let us consider the Green function for the scalar wave. 



Figure 1. Configuration of wave scattering by a black hole. 


For a monochromatic stationary wave with time dependence e the Green function 
satisfies 


- gV<l> + = - 

V~9 

where 6^{x,Xs) denotes the invariant delta function 
wave function as 


6^{x,Xs). 

- 




( 22 ) 


By separating the 


$(r,0;r„0,) = ^ 
1=0 


2 £+ 1 

Wrrs 


i>t{r.r..)P^{cos0)P^{cos0s). 


(23) 
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the radial Green function ijji{r,rs) obey^ 
(Ptjji 


drl 


+ Q(r)'0£ = -5(r* - r^*). 


(24) 


The solution of this equation can be constructed from a pair of linearly independent 
solutions of the radial equation; one is given by Min in (12) and the other has the following 
asymptotic behavior 


-^in ^ “r -^out ^ 


Wup ~ 


>'in ’ 


(25) 


r* —)■ —cx), 

e"””, r* —)■ + 00 . 

Using these two independent solutions Min, Unp, the radial Green function can be written 


as 


Mors) = 


Uin{r<) UnpM) 

W 


where r> = max(r, rs),r< = min(r, r^) and the Wronskian W is 


jjr dMup dMin . . 

W = Uin-M - Mup-l- = 2zUjAi, 

dr^ dr^ 


(26) 


(27) 


For large r, with Vg <r, the radial Green function is 

Mors) ~ _ (^_Y^2iSe _ 

2uj 

Assuming that the point source of the wave is located on —axis {6s = vr), the full 
Green function for the scattering problem is 


(28) 


<I>(r,9;r„9,) = + V)(-)’i'l(r,r,)P,(cosO). 


(29) 


£=0 


For M = 0 (free propagating wave in the flat spacetime), ipMirs) = Trsji{ojrs)M\ojr) 
and this formula reduces to the Green function of the freely propagating wave in the 
flat spacetime: 

?/ ) ^ 1 Xs| 

$0 = — ^( 2 ^ + l){Mje{(^rs)h[^\ur)Pe{cose) 


e=o 


47 r X — Xo 


(30) 


Including the term phase of the radial function (28), the asymptotic form 

of the Green function at large r is 


$(r,d) 


SivTcurr. 




£=0 




= 4>r 


giij(r3*+r*) 

dTrr, 2iu:r 


- l)Pe{cose), 


i=0 


where <Fo = *h| 5 =o represents the incident spherical wave from the point source and we 
dehned 

111 


r r r. 


f We use the relation + ^)P£{cos9)Pi{cos9s) = S{cos9 — cos9s). 


Pi{cos 6) 

( 31 ) 


( 32 ) 
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Following the standard prescription of preparing the plane wave in the black hole 
spacetime [23], the spherical wave from the point source placed at sufficiently far from 


the black hole is obtained by replacing r, r* in the phase factor of (30) to their tortoise 
coordinate r*, Vg *: 


<Fn 




2 for 9 1 and large r, Vg 


(33) 


47r(r + Tg) 

At this stage, we apply the Poisson sum formula[§]to the sum in (31) which does not 
contain the phase shift: 


J2m + i)o'-^p,{cose) 


(34) 


i=0 


+00 

2 E 

m=—oo ' 
poo 


dX Px-i/2{cos 9) e 


2im7T{\—l/2) 


H-oo 


dXXP^^Jo{Xe)+4j2 / dAAe*^Jo(A0) cos(2m7r(A-1/2)) 


m=l 


= 2iuj- 




+ 0((a;h)‘J), for 0 < 1, 

r P Tg 

where we used the asymptotic formula for the Legendre function 

B 


Px-iMcosO) 


sin 9 


Jo{X9) for A S> 1, 9 ^ vr. 


(36) 


and consider the scattering with small 9 (forward direction) in the present analysis. 
We neglect 0{u}^) contribution which is higher order in the eikonal approximation. 
Therefore, for large r, the Green function in the form of the partial wave sum is 


$ 


47r(r + Tg) 




lUTTg 


^ + M Jo((£ + 1/2) 9) 


gjaj(r,+rs*) 

AniUTTg 


^ (£+-) e*^^^e2*^Vo((£+1/2)0), for 0 < 1. (36) 


1=0 

In this formula, the factor is essential to make the inhnite sum over ^ hnite. 

Without this factor, the sum does not converge; this is well known property of the 
partial wave sum for long range potential which falls off as 1/r at large distances. We 
will discuss the issue of convergence of the partial wave sum ( [Mj ) and (36) in the next 
subsection. 


§ The Poisson sum formula is 

OO 

^nf+i/2)= E 


+00 


i^(A) 




t=o 
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We take a look at the scattering of the plane wave in the weak gravitational field of 


a point mass and check the formula (36) indeed works well. In this case, the deflection 
angle is given by the Einstein formula for light deflection 

AM AMoj 

0(4) = —- «-—, (37) 


where b ={£+ 1/2)/u is the impact parameter for the partial wave. The phase shift is 


obtained by integrating the relation (19) 




diQ{i) = —2Ma; In £ +const. 


(38) 


By taking Vg —?• oo, the incident spherical wave reduces to the plane wave. Renormalizing 


the over all factor as e'^’^”*/(Airrg) —)■ 1 and applying the Poisson sum formula to (36), 


$ = 


tur 

e 
iur 




1=0 
lujr* /*oo 


- I + 1/2)0) 


dX A X-^^^^^JoiXe) + 0{l/{ur)) 


(39) 


_ g-*2Ma;ln(4M..) ^io.r{l-ey2) - i2Mu) (2iMu, 1, 

The phase factor e-* 2 Ma;in(4Mt<;) absorbed to the constant of the phase shift. Then 

the obtained wave function has the following asymptotic form for urO"^ S> 1, 

Fh — 9iMi,A f f) \ “2+4iMa; iuj(r+2M\n{2ur)) 

* ~ iuiz-2iu:M\n(uj{r-z)) . j./r I ^ j ^ 

^ ^(l + 2^Ma;)V2; r ’ ^ ^ 

where 2 ; ~ r(l — 0^/2), r — z k, r0^/2. The first term is the incident plane wave and the 


second term is the scattering wave. Thus (39) reproduces the scattering wave function 
for 1/r potential (Coulomb scattering). For Mu: 3> l,ra;0^ = finite, corresponding to 
the scattering toward the forward direction 0 ~ 0, 

$ ~ gia7r(i-02/4) _ i2Mu) Jo{^/AMr uvO). 

2.3. Convergence of partial wave sum 


(41) 


Here we comment on convergence of the partial wave sum (36). We truncate the infinite 
sum up to f'max and the phase shift is evaluated by the WKB formula (16). As we have 


mentioned in the last paragraph in Sec. 2.1, the upper bound of the partial wave sum 
in the eikonal approximation must be £max = \/2ur, and this value becomes inhnity in 


the eikonal limit. To check the convergence of this sum, we first consider (34) which is 


5(, = Q case of (36). Let us consider the following quantity: 


-j -t-max 2 

7(9, C.x) = — J](2(' + 1) e‘‘^ J„(((> + 1/2)9). 


(42) 


f=0 


The relation (34) implies this sum becomes 


|/(9)| = l + 0((^..r)-‘) 


(43) 
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in the eikonal limit car 3> 1. We present converging behavior of this sum for 9 = 0 case 
(Fig. 2). 



Figure 2. The partial wave sum (34) as a function of fmax for ujr = 10,100,1000. 


For ~ \/2 CUV, WO can observe the absolute value of the sum oscillates around the 
true value predicted by the formula (34). This oscillation is due to interference between 
0{{ojr)^) and 0{{ur)~^) terms and the amplitude of the oscillation is expected to be 
reduced in wr 1 limit. Thus we expect that the partial wave sum with f'max = V^ojr 
converges in the eikonal limit. 

We then checked numerically the convergence of the sum of (36) for the weak 
gravitational field case that is equivalent to the Coulomb scattering. We truncate the 
inhnite sum up to the value £max and the phase shift is evaluated by the formula (38). 
The convergence of the sum (39) for 0 = 0 is shown in Fig. 3: 



Figure 3. The convergence behavior of the sum (39) for 0 = 0. 
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We expect the absolute value of this sum becomes e^^^|r(l — 2iMu)\. Fig. 3 shows the 
sum converges to the correct value if we set £max = V^uir and take cur 3> 1 limit. We 
also check the sum reproduces correct 9 dependence of the scattering wave function. 




Figure 4. Behaviors of the wave function obtained by the partial wave sum up to 

fmax(wr = 200, r = 20M). The dotted lines are the partial wave sum without the 
• (^ +1 /2) ^ 

factor e* ^ . The last panel is behavior of the analytic function that the partial 

wave sum converges in the eikonal limit. 


Fig. 4 shows behavior of the wave function (39) obtained by taking the partial wave sum 
up to iraax (tur = 200,r = 20M). For f'max = 300 ~ -s/^ur, the obtained wave function 
has a good agreement with the asymptotic analytic formula. 

For the Schwarzschild black hole case, the WKB phase shift is evaluated by the 
formula (16). 
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Fig. 5 shows appearance of the critical angular number ic ~ 51 corresponding to the 
critical impact parameter be = t'e/w = 3a/3M for null rays. To make connection with 
the geometric optics, we plot the deflection angle introduced by the relation (19): 



Figure 6. The deflection angle for null rays in the Schwarzschild spacetime calculated 
from the WKB phase shift {Moo = 10). The critical angular number is £c ~ 51. 


Near the critical angular number the deflection angle can have arbitrary negative 
values and null rays can go around the black hole arbitrary times (orbiting, Fig. 6). 
This feature is peculiar to black hole spacetimes that have event horizons. Fig. 7 shows 
the wave function calculated by taking the partial wave sum (r = 20M). 
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e 


e 


Figure 7. Behavior of the scattering wave function in the Schwarzschild spacetime 
obtained by the partial wave sum. The top left panel shows convergence of the wave 
function for 9 = 0 (r = 20M). The other panels are 9 dependence of the wave function 
with Mu! = 10, r = 20M and fmax = 100,200, 300. 


We expect the optimal value of imax for Mu = 10, r = 20M is y/2ur 280. Indeed, the 
value of the wave function for 0 = 0 is converging to the expected value for > 200 
(the top left panel in Fig. 7). 

To summarize, base on these examples investigated in this subsection, we conclude 
that the partial wave sum formula (36) is applicable to black hole spacetimes if we set 
the upper limit of the sum as £max = and take the eikonal limit car S> 1. In this 

paper, we try to evaluate this sum analytically. 


2.4- Derivation of lens equation 


Applying the Poisson sum formula again to (36), the asymptotic form of the Green 
function is 

giaj(r,+rs*) +°° roo 


$ 


Amurr. 




(44) 


Changing the integration variable to the semi-classical impact parameter b = {i+l/2)/u, 
the wave function ([44| becomes 


$ oc 



OO 




(45) 
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As the integrand is the rapidly oscillating function of b for Mu S> 1, the integral can 
be evaluated by the stationary phase method. Using the asymptotic form of the Bessel 
functionjlj] the stationary phase condition for the integrand yields 

b(- + -']=-Q(b)-2nm±9, 0(A) = 2^. (46) 

\r VsJ dX 


0(6) is the deflection function which coincides with the classical deflection angle (11). 
This is the lens equation for the black hole spacetime. For a given position of the source 
and the observer (r, r*, 0), this equation provides the impact parameter b and determines 
the trajectory of null rays connecting the source and the observer. The integer m has 
the meaning of the winding number for null rays orbiting around the black hole. The 
geometric interpretation of this equation is shown in Fig. 8: 




Figure 8. Configuration of the gravitational leasing in the geometric optics. Angles 
in the figure satisfy a + j3 = —Q + 9 in the left panel and a + (3 = —Q — 9 in the right 
panel. 


Angles a, (3 satisfies the relation 
ct + /3 = —0 ± 0, 


(47) 


where + sign corresponds to the left panel and — sign corresponds to the right panel in 
Fig. 8. Using a ^ and ^ <C 1, we obtain 


r o + r 


rr. 


= - 0 ± 0 . 


(48) 


This is the lens equation (46) with m = 0. For —0 > 27r, the null ray goes around the 
black hole more than one times (orbiting). In such a case, the winding number m is 
defined so as the condition 0 < — 0 — 27rm < 27r is satisfied. 


3. S-matrix and Regge poles for wave scattering by black holes 

For the partial waves with impact parameter b ~ 3\/3M, the scattering occurs in the 
vicinity of the peak of the effective potential (Fig. 9). 


II 


Jo (a;) 



cos(a; 


-7r/4) 


for 


|a;| > 1. 
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Figure 9. The effective potential Q{x) around the unstable circular orbit Vp. The 
quasi-normal mode satisfies Aout/^in = oo. 


Using the result of the asymptotic matching of the WKB wave function around the 
peak, the S-matrix Si = is obtained as 


A piiTU / 1 \ !^-I-1/2 

where u is determined by the relation [25] 

1 . Qo 27M^u^-e 

l' = I —; - ~ -7-. 

2 2e 

The peak of the potential for the Schwarzschild case is at Vp = 3M and 


(49) 


(50) 


Q« = Q(rf), Q; = ■ (51) 

Absolute value of Si represents the reflection rate of incident waves. As the black 
hole absorbs waves, the phase shift acquires the imaginary part and the reflection rate 
becomes smaller than unity for small impact parameters. Indeed, 




= e 


— 2(5j 


1 


1 + exp 




1 

1 + exp (—vr(£ — ^c)) ’ 


(52) 


where ^c = 3\/3Mu is the critical angular momentum. The waves with i < ic are 
absorbed by the black hole, be = = 3\/3M is the critical impact parameter 

(Fig. 10). 



Figure 10. The reflection rate for = 10. 
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The S-matrix (49) has poles in the complex i plane of which location is determined 


by poles of the Gamma function r(— z/), z/ = n, rz, = 0,1, 2, • • That is, poles of the 
S-matrix in the complex i plane are 


C ~ 4 + * , 4 = sVsMu, 


rz, = 0,1,2, ■ ■ ■ 


These are Regge poles 


In the vicinity of the rzth pole £n, the S-matrix is 


S, 


In . / 


In = (-) 


n+1/2 


,-(n+l/2)_ 


£-£n ^ ^ 4^ V4y V” ' 2y n\ 

For rr 3> 1, rz! Ri y/27in e“”rz"', (1 -|- l/(2rz))"'+^/^ hence the residue 7 „ is 


In 


^ ’ 2ttL' 


(53) 


(54) 


(55) 


4. Evaluation of scattering wave 


The formula of the Green function (44) is decomposed to rrz = 0 and rrz 4 0 parts: 

git*;(r,+rs*) 


$ = 


Airiurr. 


Jo(A0) + 




d\ A e^W-i/2 j^(^xe) 


As we have discussed in Sec. 2.3, the integer m is the winding number for null rays 
orbiting around the black hole. Hence the hrst integral term represents contribution of 
rays which directly reach the observer (direct part). By identifying the phase shift as the 
lens potential of the gravitational lensing, the direct part reproduces the Kirchoff-Fresnel 
diffraction formula for the gravitational lensing HD]. On the other hand, the second 
integral term represents the contribution of rays go around the black hole (winding 
part). We evaluate these two parts separately. 


.( 66 ) 


y.i. Winding part : contribution of Regge poles 


As the S-matrix has poles in the complex i plane, the integral of the second term in (56) 


can be evaluated by applying Gauchy’s theorem. We use the contour shown in Fig. 11. 
We hrst show that the integral along the circle |A| = i? vanishes as i? —?• cx). Introducing 
a new integration variable = —'zA^/(2co’r) 


dAAe*^e2*'^"-i/Vo(A0)e 


227rm(A—1/2) 


J\\\=R 

^ The quasi-normal frequency is given by treating £ as real: 

1 


343M 343M r ' 2 


To show that the integral along the circle vanishes, we use the formula 

9 1 r+°° 9 

e-y=^ dke-’^+^^^y. 

\IJ — f-v -1 
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'\X\=R 


dX \exp 


2ke 


, _+ i27rm ) A 
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where we choose a = 57r/4 for Im(A) > 0 and a = —vr/d for Im(A) < 0 for convergence 
of the integral. Hence, as i? —>■ oo, the integral becomes zero. 



Figure 11. Poles of the S-matrix are situated in the first quadrant of the complex A 
plane. 


Thns the second integral in (56) has contribntions from poles and the integral along the 
imaginary axis. 

oo cxo 2 ^—iOG 

EE 

m=l n=0 

The integral along the positive imaginary axis is 


JO 


m=l 


m=—l ' 


~ p+ioc p+ioo ■ i7r(A-l/2) o 

^,Jo k 2sm(7r(A-l/2)) ' 


m=l ' 


—ttA 


= 0(0;'^), 





„e jQ(iA0) 


(58) 


and yields higher order contribution in the eikonal limit because the hrst integral in (56) 


gives 0{u) contribution. The contribution of the integral along the negative imaginary 
axis is also 0{u^). We ignore these terms in our treatment. Thus, the wave function 


(56) becomes 


4> 


gjaj(r*+rs*) 

Airiurr. 


dX X Jo(A0) + 27rt 4 7n MLO) /(4) 


n=0 


where 


/(4) = E 


227rm(£ri —1/2) _ 


^2i7Ti-n 


( 60 ) 
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The sum over n can be evaluated by applying the Euler-Maclaurin formula and keeping 
the contribution of the leading order in the eikonal limit; 

CO 2 

2z7r ^ In MLO) /(4) 


n=0 


Airlc 


dn { £c + - + in ] e 


§+*”■) J, 


ic + - + in ] 6 


V 2 Vvcur 


We have used 


/(^=+2 


~ _e2i7r£cg-7r^ 


After all, we obtain 

giaj(r,+rs*) 


$ 


Aniurr 


s L^o 


/ I + 2 + “ 


(61) 


d\\e‘^r e“^> Jo(A9) + ,/- £-”-'■’■'■‘+“’‘'63* 


(62) 


(63) 


4-2. Direct part and total wave function 

We evaluate this term assuming that impact parameters of corresponding null rays are 
not so small and gravitational field can be approximated by the Newtonian weak field. 
Within such an assumption, the phase shift can be approximated by that of the weak 


field form (38) and the direct part can be obtained as (39). Therefore, the formula for 


the total scattering wave is summarized as follows: 

gitj(r*+rs,) p 


$ 


47r(r + Vg) 


ciJo{bEUj 6 ) + C 2 Jo{bcUj 6 ) , for 0 < 1, 


(64) 


Cl = e 


—i2Mu In 2Muj ^nMu 


e^^W{l - i2Moj)e—^ 


-4ajf6l2 


C2 = 


TT 

2 ljjf 




where = \/AMf is the Einstein radius and be = 3a/ 3M is the critical impact parameter 
for null rays. This formula is applicable for r, r* S> 2M and Mu S> 1. Let us assume 
that the point source is near the black hole and the observer is sufficiently far from the 
black hole (r ^ r^). Then r ^ Vg. Of course the location of the source must be far 


from the black hole to apply (64), we dare to consider such a situation to investigate 


the wave optical image and the interference effect by the black hole. As the impact 
parameters for the direct rays must be larger than those of the winding rays, we must 
require 6^ > be which constrains the distance of the point source from the black hole in 
our approximation: 


27 

Vg > —M = Tc ~ 6.8M. 

s 4 c 

The ratio of two coefficients is 

|C2| 3 im 

-\l -X e = e 

Cl 2 


- 1/2 


0.043 X 


6.8M 


- 1/2 


(65) 


( 66 ) 
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In the geometrical optics, it is known that the intensity of the winding ray reduces to 
the factor ~ 0.0019 compared to the direct ray [5]. In terms of the amplitude of 
the wave, this value corresponds to e~'^ ~ 0.043 and is consistent with the obtained 
value |c 2 |/|ci| in the eikonal limit of the scattering wave function. 


5. Applications 


As applications of the analytic formula (64) of the scattering wave, let us consider images 


of the black hole and the interference effect in the frequency domain. 


5.1. Wave optical image of black holes 


Using the scattering wave function (64), we can reconstruct wave optical images of black 
holes illuminated by a point source^ . As the imaging system, we consider a detector 
with a convex lens which transforms the interference fringe to the image (Fig. 12). 



Figure 12. The configuration of the imaging system with a convex lens. 


The location of the lens of the detector is Z = 0 in the detector’s coordinates and the 
wave function just in front of the lens is ^(X, Y, 0). Using the Fresnel-Kirchoff diffraction 
formula [26], the wave on the detector’s focal plane Z = f is given by 

$7(X/,y/)oc/ X ^ (67) 

Jx^+Y^<a2 ^ 


where X/, Yj are coordinates on the focal plane, a is a radius of the lens and £ is the path 
length between a point on the lens plane and a point on the focal plane. The function 

X^ I y ‘^ 

2 / in (67) represents the action of the convex lens which transform incident plane 
waves to spherical waves focusing at the focal point. Assuming that the focal length is 
sufficiently larger than the radius of the lens. 


= V(^ - XiY + (u - y»2 + p^f + 


(X - Xjf + {Y- Yj 

2 / 


( 68 ) 


and we obtain 


^i{Xj,Yj) oc [ dXdY^{X,Y,0)e~^7^^^^+^^^\ (69) 

J X^+Y^<a2 
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Thus the wave function on the focal plane is the Fourier transformation of the wave 
function <F. The relation between the world coordinates {x, y, z) and the detector’s 
coordinates (X, Y, Z) is 

X = r sin 9 cos 0 = X cos 9q — Z sin 6*o, 

y = r sin 6 sin (j) = Y, (70) 

z = r cos 6 * = X sin 6*0 + 2’ cos 6 *o, X^ + Y‘^ + Z'^ = r^. 


It is possible to obtain cos 6 * as a function of X, X: 


cos 6 * 



sin 6*0 + 





cos 6 * 0 - 


(71) 


Fig. 13 shows an example of images of the Schwarzschild black hole illuminated 
by a point source. Images are obtained by drawing |<h 7 -| which is the amplitude of the 
Fourier transformation of $. 



Figure 13. Images of the black hole for different viewing angles Oq = 
0.0, 0.02, 0.05,0.09 (from the left to the right panels). Parameters are = lOM, r = oo. 
Moj = 800 and a = 0.03. 


In our previous paper [9], we numerically solved the Helmholtz equation in the 
Schwarzschild spacetime and images of the black hole are reconstructed using the 
obtained wave function. Due to limited resolution of the numerical calculation, it was 
difficult to perform simulation with high frequency waves corresponding to the eikonal 
limit and it was hard to resolve the fine structure of black hole images. In the present 
analysis, as we have obtained the analytic form of the scattering waves, it is possible to 
acquire images using high frequency waves. In Fig. 13, we can observe a bright outer 
ring which corresponds to the Einstein ring by direct rays and a thin inner ring which 
corresponds to the rim of the black hole shadow (photon ring). It was not easy to 
recognize this structure of the black hole images in our previous paper [U]. 


5.2. Interference effect in wave scattering by black holes 


We can examine the interference effect in the frequency domain using (64); that is caused 
by the interference between the direct rays and the winding rays. For bEUj9,bcUj9 ^ 1, 
using the asymptotic form of the Bessel function, the scattering wave (64) can be written 


as 


$ ~ di cos(6^a;6*o) + ^2 cos(6ci^6*o). 


( 72 ) 
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where di, d 2 represent numerical factors independent of u with their ratio 


Mil 



(73) 


After averaging the amplitude of the wave over rapidly oscillating frequency scale 
~ l/{M6o), the power spectrum is 


I{oj) = |<h| 


MiP , M: 


+ 


P did*2 + d\d2 u \ Q 

- 1 - ^ -cos [[pE — OcjOJOQ 


(74) 


The first term and the second term represent the intensity of the direct wave and the 
winding wave, respectively. The third term represents the interference between the direct 
wave and the winding wave. For the point source at ~ 7M, the impact parameters 
for the direct rays and the winding rays become nearly equal and the interference term 
results in the “beat” in the frequency domain. The modulation due to the interference 
appears in the power spectrum (Fig. 14). 



Figure 14. The power spectrum of the scattering wave with = 7M. The amplitude 
of the modulation is A/// ^ 0.1. The period of the modulation depends on the viewing 
angle 6q. 


The period of the modulation in the power spectrum is 

-1 


MAu = 


271 1 


do sVS VV 27M/4 


- 1 


(75) 


As the typical value for the viewing angle ^o, if the Einstein angle 6 e = \/~^Mjr~s is 
assumed, we obtain 


MAu 


TT 


2-3 



r., > 


27M 


(76) 


For Vs ~ Tc, the visibility of the modulation is 

y ^ /max - /min ^ ^1^2 + ^ 1^2 _ ^ _ q 

■ /max + /min W + W 

This value is not so small and we may have the possibility to detect this interference 
effect observationally, which means direct verification of the black hole spacetime. 




galactic core BH 

intermediate mass BH 

stellar mass BH 

mass 

IWMq 

lO^Mo 

lOMo 

Aw 

40Hz 

40kHz 

4MHz 


Table 1. The typical frequency of the modulation in the power spectrum. 
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Table I shows the frequencies Au for black holes with different masses. For intermediate 
mass black holes and stellar mass black holes, the modulation in the power spectrum is 
expected to be detectable using sub-mm radio telescopes with suitable band width and 
the frequency resolution. 

6. Summary 

We investigated the wave optics in the Schwarzschild spacetime following the standard 
treatment of the wave scattering problem and evaluated the Green function for the 
stationary monochromatic point source in the eikonal limit. To make the partial wave 
sum of the scattering wave converge, we retain the next to leading order contribution in 
the phase of the partial waves for large r. Effect of the orbiting scattering is taken into 
account as the contribution of Regge poles to the scattering wave. The wave optical 
images of the black hole are obtained and the interference effect in the frequency domain 
is discussed. 

As the straightforward extension of the analysis performed in this paper, the 
formulation of the wave optics in the Kerr spacetime is now on going. We expect 
that the effect associated with the black hole spin such as the superradiance adds new 
features to the scattered wave. We will report on this subject in our forthcoming paper. 
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